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ABSTRACT 



We present here the solution for accretion disk structure, which describes 
continuously the transition between optically thick and optically thin disk 
regions. We show that the disk structure equations without advection, used here, 
give two branches of solutions which do not intersect for L < Lb <^ 0.6ivEdd for 
a = 1.0 and Mbh = 10 8 Af e . For larger luminosities there are no global solutions 
of the equations without advection. We suggest, as one of the possibilities, that 
advection becomes important even before, at luminosities ~ O.livEdd, when the 
inner disk regions become unstable due to radiation pressure dominance and 
induce the transition from optically thick to optically thin disk in the global 
solution. 



Subject headings: accretion, accretion disks - black hole physics - radiative 
mechanisms 
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1. Introduction 

It is well known that disk accretion plays a crucial role in many astrophysical objects 
and phenomena. Linden-Bell (1969) was the first who proposed the model of gaseous 
disk accretion onto black holes as a source of energy of quasars and active galactic nuclei. 
Shakura (1972), Pringle & Rees (1972) and Shakura & Sunyaev (1973) built the Newtonian 
models of the accretion disk, while Novikov & Thorne (1973) and Page & Thorne (1974), 
gave the theory of disk accretion in the framework of General Relativity. This theory is 
important for the structure of the innermost parts of the disk. 

In the subsequent period a lot of work was devoted to this problem. In the most early 
papers (see e.g. Lamb 1991 for a review), the emphasis was mainly on the existence and 
analysis of the local disk solutions at a given radius, r. Of greater importance, however, is 
whether or not these local solutions at different radii connect smoothly to form physically 
self-consistent global disk solution that is continuous from large radii to the inner disk edge. 

It is known that for given values of the accretion rate, M , and viscosity parameter, 
a, there are two different solutions of the steady state disk structure equations at a 
given radius, r. One of them corresponds to an optically thick, low-temperature state 
(T ~ 10 5 — 10 7 K; Shakura & Sunyaev, 1973). The other solution corresponds to an optically 
thin, high-temperature state (T ~ 10 9 K; Shapiro, Lightman & Eardley, 1976). 

The physics and properties of these solutions, e.g. stability, are quite different. The 
formulae for radiation flux and radiation pressure conventionally used are valid only in 
limiting cases either for the optically thick region, or for the optically thin one, but not 
for the intermediate zone. So far, only a few papers discussing the intermediate zone have 
been published (Maraschi & Molendi 1990, Hubeny 1990, Liang & Wandel 1991, Wandel & 
Liang 1991, Lasota & Pelat 1991, Kusunose & Mineshige 1992, Luo & Liang 1994). 
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Here we study the disk behavior at all r, including r e jj 1, and derive a new 
approximate formula for the radiative flux which is valid at all optical depths. 

Using this formula, we follow Bjornsson & Svensson (1991) for determination of the 
equilibrium structure of the geometrically thin disk and analyze the general topology of 
the family of the global disk structure solutions for a fixed viscosity parameter a. Each 
curve corresponds to a given rate of accretion, M. From this analysis, a bifurcation of the 
solutions into different sub-families is clearly seen. In order to isolate the main behavior of 
the global solutions we analyze the simplest accretion disk model, where the radiation field 
in the optically thin limit is pure electron-proton bremsstrahlung. 

In Section 2 we derive an approximate formula for the radiative flux valid for all r and 
in Section 3 we present the disk structure equations. In Section 4 we discuss numerical 
solutions of the disk structure equations. In Section 5 we discuss the applicability of our 
results to realistic disk models and give a qualitative descriptions of the realistic accretion 
disk structure at high luminosities, where advection may be critically important. 



To obtain the formula for the vertical radiative flux, -F ra a, in the disk and the radiative 
pressure, P ra a, for any optical depth, we consider the Eddington approximation in a grey, 
scattering medium in LTE conditions. The first momentum of the equation of transfer for 
plane-parallel geometry gives 



where z is a vertical coordinate (z = 0, at the symmetry plane of the disk), p is the matter 
density, S is the energy density and Q is the source function. Using the mean value of p } 



2. 



General Formula for the Radiative Flux 





- 5 - 



one can write down the following approximate formula for the radiation flux 

i^rad = (2) 

where F is the surface value of -F ra a, F =F TeL d\ z =h } h is the half height of the disk and 
Yj =2ph } is the surface density. Using equation (f) in equation (2) we get the expressions 
for S and P ra a in the Eddington approximation, 

3P - = s ' = " r4 ( 1 -^)- < 3 > 

From the second momentum of the transfer equation we obtain 

c — - — = -npF rad = -2—pnpz, (4) 
az L 

where, n = crj/m p} aj is the Thomson scattering depth, m p is the proton mass and c is the 
speed of light. Equation (4) can be rewritten in the following form (using r = npdz } p 
is assumed to be constant, r = r — npz, where r = uph). 

= F rad = 2^(to -r), (5) 
(It kIjq 

that we solve with the boundary condition: 

j-, I _ 7? _ C S\t=0 _ 3cP ra d| T =0 

-Trad|r=0 — f — ~ — ~ , 

resulting in 

cPradlro = ~ T~ = + ~ ^T" = F 0{^ + Tr), (6) 

where P ra( i| To is P ra d in the mid-plane of the disk, where r = r . 

Eliminating S from eq's (3) and (6), we obtain an expression for the temperature in 
the mid-plane of the disk T c = T\ T=0 

T 4 3 F f 4 2 \ 
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Here we have introduced the total optical depth to absorption, r a <C r, 

aT c 4 c 

Introducing the effective optical depth 

we find from equation (7) the connection between the flux from the disk surface, F 0} and T c 

2aT*c / 4 2 \ _1 

— — H 1 

3r V 3r 3t2 



^0 = ^(1 + ^ + ^ • (9) 



Equation (9) has the following properties: 

1) In the optically thick case, r 3> 1, r* ^> 1, and we have the diffusion limit 

P 2aT c 4 c 
3r 

2) In the effectively optically thin case, t* <C 1, any r 3> r*, we have 

= Qph. 

Thus equation (9) can be used for any optical depth. 

Our consideration is in part similar to the one of Hubeney (1990), but has more general 
character, because we have taken electron scattering into account, which is the main input 
in the pressure balance equation (eq. [4]) and gives negligible input in the energy equation 
(eq. [1]). In our case n » kjj, whereas Hubeney (1990) considered only absorption with 
k = Kff. In addition we used the Eddington boundary condition, while Hubeney used 
Chandrasekhar's ones. 

The main differences, however, is that we use the a-description of the viscocity, and 
close the system by our eq. (2) (energy production is proportional to the mass). Both 
approximations have been introduced by Shakura (1972) and are of common use. On the 
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other hand, Hubeney (1990) used the viscosity approximation of Lynden-Bell and Pringle 
(1974) and an isothermal approximation for the z-disc structure. 

His final formulae therefore are quite different and are much more complicated. Our 
approach seems to be more physically consistent and more general. 



3. The Disk Structure Equations 



With equation (9) we can investigate the geometrically thin disk structure without 
any restriction on the optical depth. To that equation we add the three disk structure 
equations: (i) Hydrostatic equilibrium in the vertical direction 

Ptot = —n p (-) r;\ (10) 
r g a T \rj 

where, r* = r/r a , is the dimension-less radial coordinate, r g = GM/c 2 } and where we have 
introduced the dimension-less proton concentration, n p = n p r g crj. (ii) Conservation of 
angular momentum is expressed by 

aP tot = ^m(-) \-JI 2 J(r^ (11) 

where, J(r*) = 1 — (6/r*) 1 / 2 , accounts for the stress free boundary condition at the inner 
edge of the disk, and a is the viscosity parameter of Shakura (1972), Shakura & Sunyaev 
(1973). The dimensionless accretion rate is, m = Afc 2 /ivEdd, where i^Edd is the Eddington 
luminosity, (iii) The rate of dissipation of gravitational energy equals the local radiative 
cooling rate that in turn escapes as local vertical flux, providing us with the local energy 
balance of the disk, 

E+ = £ rad = (y) , (12) 



where the dissipation rate is given by 



3m„c 3 . / h x 1 



E + = —?-m - r; 4 J(r,), (13) 



- 8 - 



while the radiative cooling rate can be expressed as 

3a f m e c 3 2 , N 

^rad = ' n\F hr . 14 
87rr^(7T 

Here, ctj = 1/137, is the fine structure constant and Ff, r is the dimension-less cooling rate 
due to bremsstrahlung. In our numerical calculations below, we restrict Ff, r to include 
electron-proton bremsstrahlung only. 

Our equation of state contains contribution both from gas pressure as well as radiation 
pressure 

-Ptot = -Pgas + -Prad, (15) 

where the gas pressure is given by 

„ m P c 2 ^ 

P gas = 2^n p 0, 16 

where m e is the mass of the electron and the dimension-less temperatures is defined as 
= kT/m e c 2 . The factor 2 arises as we consider contributions to the gas pressure from 
both the protons and the electrons, that in this work are assumed to be at the same 
temperature. 



The expression for the radiation pressure, with account of equations (6) and (9), i 



is 



„2 \ 4 



rad 3 \ k " J ( ^ 3r ) ( ^ 3r ^ 3t t„) ' ^ ^ 

and the free-free optical depth is evaluated according to (using eq. (8) and the right hand 
equality of eq. (12)) 

To = Ijl(J^) (^-)r*E rad e;\ (18) 
ac \m e c z ) \r J 

As the expressions for the flux, F 0} and the radiation pressure, P ra a, are rather 
complicated functions of optical depths and temperature, an analytic disk solution for the 
general case is not possible. 
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4. Results 

We solved the system of equations discussed above numerically, using the technique 
described in Bjornsson and Svensson (1991). We introduce a local compactness parameter 
based on eq. (9) that bridges the optically thick and thin conditions. We then solve 
equations (7) and (8) in Bjornsson and Svensson (1991), modified to include other details 
of our models. 

In Figure la we plot the optical depth, function of radius, r*, for the case 

Mbh = 1O 8 M , a = 1.0 and for various values of ra. In Figure lb we plot the analogous 
dependence of the dimension-less temperature, 0, on r*. One can clearly see the bifurcation 
of the global solutions into optically thick and optically thin families. For low accretion 
rates ra <J 10~ 2 our numerical solutions agree very well with the analytical solution of 
Shakura & Sunyaev (1973) in the optically thick case, and Bjornsson & Svensson (1992) in 
the optically thin case. 

For accretion rates ra <J 9.4, there are thus two families of the global solutions which 
are continuous from large radii to the inner disk edge. At ra 9.4, the two global solutions 
cross each other around r* = 32/3, and for this case we observe the transmutation of the 
optically thick solutions into the optically thin one and vice versa. For ra ^ 9.4, there are 
no global disk solutions. 

For each fixed ra ^ 9.4 there are two formal branches of solutions which both start and 
end either at large radii or at the inner disk edge but without any solutions for r* in some 
region around r* = 32/3. At the turning points of these formal solutions, the gradient of 
all physical quantities (as functions of r*) tend to infinity, violating the assumptions of the 
geometrically thin disk equations. 

The value of 9.4, is the maximum possible value of the rate of accretion for 
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which global disk structure solutions exist. Note that the bifurcation point is the point of 
maximum energy release in the disk. It is located at the radius r* = 32/3 that maximizes 
r* 3/,2 J(r*) and is independent of the viscosity parameter a (see e.g. Bjornsson & Svensson 
(1991) for a discussion). In Figure 2 we plot the dependence of mj on the value of the 
viscosity parameter a. 

It is obvious that mj can be both greater and smaller than the value of ra corresponding 
to the Eddington luminosity. Our calculations verify also that mj is very insensitive to the 
mass of the black hole. Let us stress that our results reflect the formal properties of the 
accretion disk equations in Section 3, which have physical meaning only at m < 1/e for a 
global solution to exist, or for sufficiently large values of r* in the case m > 1/e. Here, e is 
the efficiency of the accretion, e = L/Mc 2 } e ~ 1/16. Formally, at m > 1/e, the ratio h/r is 
becoming large in the inner parts of the disk and the disk approximation fails. The increase 
in disk thickness indicates the beginning of intensive mass loss from the inner parts of the 
disk in combination with the disk accretion further out (Bisnovatyi-Kogan & Blinnikov 
1977). 

In a subsequent paper we shall discuss extensions of this work, including more detailed 
micro-physics like the Comptonization effect, two-temperature structure, electron-positron 
pair production and advection. 

5. Discussion 

For making a choice between the optically thick and optically thin solutions obtained 
above (in order to apply them to accretion disks existing in reality) we must use stability 
criteria. As was shown (see e.g. Pringle, Rees & Pacholczyk, 1973; Lightman & Eardley 
1974), the regions of the optically thick disk dominated by radiation pressure are thermally 
and secularly unstable. The onset of the thermal instability is seen in Figure 3, where 
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the optically thick curve makes a transition from the positive slope (corresponding to the 
stable gas pressure dominated state) to the negative slope (corresponding to the radiation 
pressure dominated unstable state). 

A consequence of the instability development is generally believed to be the transition 
from the optically thick (geometrically thin) to the optically thin (geometrically thick) disk. 
This transition is accompanied by the buildup of a large pressure gradient and angular 
velocity drop that make it necessary to include advection in the equations for describing 
this transition region structure (Bisnovatyi-Kogan & Paczyhski, 1981). We suggest that 
in presence of advection (of entropy in the energy equation and velocity in the equation 
of the radial motion) a third branch appears, in addition to the two separate branches of 
solutions, presented here (Figure 4). The third branch describes the transition from the 
optically thick (outer region of the disk) to the optically thin (inner region of the disk). This 
transition is shown schematically by the dash-dot curve in Figure 4. The rise in pressure 
during this transition must be accompanied by a corresponding drop in angular velocity, so 
the resulting acceleration in the transition region is directed towards the black hole. 

The onset of the instability and creation of a disk with optical depth transition happens 
when the luminosity of the accretion disk is about 0.1 of the Eddington critical luminosity 
(Shakura & Sunyaev, 1973; Bisnovatyi-Kogan & Blinnikov, 1977). This happens when the 
radiation pressure dominated region is considerably developed in order to transform the 
local instability, discussed above, into a global one. Subsequent increase of the accretion 
rate must not change qualitatively the structure of the disk. 

Our solution, where the advective terms are not included, has a critical luminosity 
L = Lj, 0.6ivEdd (which corresponds to mj 9.4), above which the global continuous 
solutions cease to exist. 1 



1 Note that on this luminosity the matter outflow begins from the accretion disk caused 
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We believe that the existence ol the bilurcation point is a direct result of neglecting the 
advection terms. A more carefully calculated solution (including advection) with the optical 
depth transition, should not feel this point and should smoothly extend into the region 
L > Lb. The non-physical solution with the infinite gradients with respect to r, should not 
exist in a more detailed treatment. The behavior of the solutions neglecting advection, and 
their suggested structure including advection are indicated schematically in the Fig. 4. 

We want to emphasize that for accretion rates as high as we consider in our analysis, 
the conditions in the optically thin plasma requires more complex physics to be included 
(e.g. two temperature plasma, comptonization, pair creation). The goal of our analysis, 
however, was to emphasize the principal properties of the solutions of the structure of the 
disc as a whole, which do not depend on physical detailes. 

Finally we would like to note, that the qualitative conclusions discussed above need to 
be checked by solving the equations including advection and a continuous description of the 
optical depth transition (see formula (17)), a work currently in progress. 

We acknowledge prof. M. Abramowicz for useful discussion. This paper was supported 
in part by the Danish Natural Science Research Council through grant 11-9640-1, in 
part by Danmarks Grundforskningsfond through its support for an establishment of the 
Theoretical Astrophysical Center. G.S.B.-K. acknowledges the partial support of RFFI 
grant 93-02-17106, and Astronomy Programm 3-169. 



by its radiation field (Bisnovatyi-Kogan & Blinnikov, 1977 ). 
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FIGURE CAPTIONS 



Fig. 1. — a) The dependences of the optical depth r on the radius, r*, for the case 
Mbh = 10 8 Af e , a = 1.0 and various values of m. The curves are numbered in order of 
increasing m = 1.0,3.0,8.0,9.35,10.0,11.0 and 15.0. The upper curves correspond to the 
optically thick family, the lower curves correspond to the optically thin family, b) The 
temperature structure for the same model as in a). The curves are numbered as in a). 

Fig. 2. — The dependence of the maximum attainable accretion rate, mj, on the value of 
the viscosity parameter a. 

Fig. 3. — The accretion rate, Tfi , clS 9j function of the Thomson scattering depth r for the 
same model parameters as in Fig. 1, evaluated at r* = 32/3. The curves are numbered as 
follows: 1 - the optically thin branch, 2 - the optically thick radiation pressure dominated 
branch, 3 - the optically thick gas pressure dominated branch. 

Fig. 4.— The total pressure Ptot &s ci function of the radius r* for a = 1.0. The solid curve 
corresponds to the accretion rate m = 9.35, the dotted curve corresponds to m = 10.0, 
and the dash-dot curve shows schematically the solution, taking into account the effects of 
advection (see text for further discussion). 



